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ABSTRACT 


A new algorithm for calculating h” optimal controllers 
is investigated. The new algorithm is significantly simpler 
than existing approaches and yields much simpler 
controllers. The design equations are first presented. 
Special system transformations required to apply the new 
algorithm are then presented. The use of the new algorithm 
with sampled-data systems is outlined in detail. 

Several constraints on the characteristics of the 
problem formulation are required for the application of the 
design equations. The consequences of these constraints are 
investigated by applying the algorithm to a simplified 
design for a subsystem of a large space structure ground 
test facility. The investigation of these constraints is 
continued by application of the design equations and 
constraints to an extremely simple tracking problem. The 
result of these investigations is the development of a 
frequency dependent weighting strategy that allows realistic 
control problems to be cast in a form compatible with the 
new algorithm. 

Further work is indicated in the area of developing 
strategies for choosing frequency-dependent weights to 
achieve specific design goals. The use of the freedom in 
problem formulation to achieve robustness/performance 
tradeoffs should also be investigated. 

It is not clear that the new algorithm always leads to 
simpler controllers. The more restrictive formulation may 
dictate that frequency-dependent weighting adds to the 
controller order disproportionately. This effect must also 
be investigated. 
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NOMENCLATURE 


<W A > 


^x(A) 


F ( s ) 


IMC 

BET 

LOS 

LQG 


minimum singular value of a matrix A 

maximum singular value of a matrix A 

infinity norm of a system transfer function 
matrix 

image motion compensation 
base excitation table 
line of sight 
linear-quadratic Gaussian 
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INTRODUCTION 


Until the recent work by Glover and Doyle m , the design 
of H w controllers promised to be a long and arduous ordeal 
for the designer. Moreover, the resulting controllers 
tended to be extremely complex, sometimes exceeding the 
order of the control model by a factor of five. Their 
publication of design equations for controllers of the same 
order as the control model is thus a significant advance in 
the state of the art. However, the question of whether H* 
control techniques can be successfully applied to large 
space structure (LSS) control design problems is by no means 
answered . 

The foremost question in the mind of any LSS control 
designer is that of applicability of h" techniques to the 
usual goals of LSS control. It is known that H* 
optimization can achieve (at least mathematically) any of 
the goals of disturbance rejection, command tracking, and 
robustness. What is not known is whether h” design methods 
can be used to design controllers which simultaneously give 
acceptable performance and do not suffer from the known 
shortcomings of LQG techniques, e.g., lack of robustness. 

The purpose of the work presented here is to address the 
issues of the applicability of the new algorithm for 
calculating h“ controllers for large space structures. 

The organization of the report is as follows. Section 
1.0 contains a brief discussion of the h“ performance 
criterion and the design equations which must be solved in 
order to find the optimal H° controller. 

Section 2.0 describes a transformation required to 
satisfy constraints on using the design equations. The 
result is a set of equations that can be readily used to 
transform a given state space realization to one of the 
required form. 

Section 3.0 outlines modifications to the state space 
formulation that are required in order to apply H design 
formulas to sampled data systems. The result is a set of 
state space formulas for applying the well known w-plane 
transformation to multivariable control problems. The 
equations, minus the derivation, can be found in Glover. An 
outline of the derivation is included here for completeness. 
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Section 4.0 documents the application of H* techniques 
to a simplified model of the ACES IMC subsystem, including 
BET disturbance effects. Section 5.0 uses a simple tracking 
problem to discuss weighting schemes which allow the new 
algorithm to be applied to realistic control design 
problems . 

Section 6.0 contains conclusions and recommendations 
for further work. In particular, it is suggested that 
strategies for developing frequency-dependent weights be 
investigated . 


1.0 h" performance criterion and design equations 

The h” control problem can be stated as follows. The 
system equations are 


x = Ax + B,w + B 2 u 


z = C.,x + D n w + D 12 u 


y = C 2 x + D 21 w + D 22 u 

where w is in R m1 , u is in R 1 " 2 , z is in R p1 , and y is in R^. 
The signal w is an exogenous input which may be either 
disturbances or command signals; u is the control input 
vector; z is actually the performance related vector; and y 
is the vector of measurements that is actually available for 
feedback. 

Although most calculations with this technique are done 
in state space form, the performance criterion is most 
easily stated in the frequency domain in terms of the closed 
loop transfer function matrix. The open loop transfer 
function matrix can be expressed in terms of appropriate 
partitions as 


[A, B, , C,, Djj ] 
g , 2 (s) = [A, B 2 , C,, D 12 ] 
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^21 ( s ) — [A, t ^2' ^21 ^ 


G ->2 ( S ) [ A , Bj r ^2 1 ^22 ^ * 


If a controller with transfer function matrix K(s) is 
connected from y to u, the closed loop transfer function 
matrix is given by 

T(s) = G„ + G 12 K(I - G 22 K)' 1 G 21 . 


The h” control problem is to find a controller K which 
yields a stable closed loop system and for a prespecified 
real number r, 


I! t | | ffl < r 


where 


I T | L = S w P a max (T ( j w) ) • 


In the special case of a scalar transfer function, the goal 
can be stated as simply insuring that the closed loop 
frequency response has a magnitude less than r. This is 
clearly important in disturbance rejection problems and can 
be used also for tracking type problems. 


The H° optimal control problem is that of finding the 
smallest such r such that a stabilizing controller exists. 
Note that once a method of satisfying a given T bound is 
identified, the job of obtaining an optimal solution is not 
difficult, although it is iterative. The significance of 
the recent work of Glover and Doyle is that their equations 
yield not only a controller that achieves the prespecified 
bound but is only the order of the original control model. 
The original factorization algorithms tended to yield 
controllers many times as large as the original plant. 

Several constraints must be placed on the plant 
equations in order to apply the design equations. The first 
constraint is that the realization (A, B 2 , C 2 ) be 
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stabilizable and detectable, as is usual. 

The second constraint is actually two constraints that 
are artifacts of the derivation procedure and are required 
for well-posedness . They are 


rank D 12 = m 2 


and 


rank D 21 = p 2 . 


One of the consequences of the above rank conditions is that 
the control input u must appear in the performance oriented 
regulated variable z. This condition is similar to a 
condition required for well posedness of the LQG problem. 
These conditions also insure the calculation of a realizable 
controller. It will also be seen in subsequent sections 
that these requirements cause some difficulty in calculating 
a low order controller. 

An assumption which is not independent of the above 
rank conditions is that 


d 12 = [0 I] T 
and 


D 21 - [0 I ] . 


Transformations to achieve the required forms for D 12 and D 21 
will be derived later. 

The final constraints are necessary for the design 
equations to yield a solution. They are not necessary in 
the strict sense. However, it is not known how to calculate 
controllers to achieve the bound when these conditions are 
violated. Sufficient conditions for the design equations to 
work are 
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The realizations of G 12 and G 21 are minimal. 

rank G 12 (jw) = m 2 for all w 

and 

rank G 21 (jw) = p 2 for all w. 

The design equations can be presented in an abbreviated 
form by denoting the solution to the Riccati equation 

Q + XA + A T X - XPX = 0 

by its Hamiltonian matrix 


X = Ric 


A -P 


-Q -A' 


D n is also partitioned as 




D nn 

D 1112 

D n 


D 1121 

D 1122 

Two 

intermediate 


variables 


are defined: 


D 1x t D 11 D 12^ 

- [D T n d t 12 ] t . 

Then define 
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R 


= D U D U 



0 


0 

0 


D xiD T x1 


r %i 0 

0 0 


X m = Ric 


U -P 

T ! 

-Q -U T 


where 


-1r>T 


P = BR 'B', U = A - BR" 1 D t 1 x C 1 , and Q = - C l ,D 1 x R' 1 D l 1 x C 1 


-i^t 


Similarly, 


Y„ = Ric 


U -P 
-Q -U T 


where now 


P - C^C, U = A T - C t R‘ 1 D x 1 B t 1 , and Q = B 1 B T 1 - B^^R'^B^ . 
Similarly to LQG design, matrices F and H are defined as 
F = -R 1 [ D t 1 x C 1 + B t XJ = F \ 2 F T 2 ] 
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H = -[B,D T x1 + Y ffl C T ] R 1 = [H n H 12 H 2 ]. 

At this point, it is necessary to perform tests to 
determine whether it is possible to achieve the specified r 
by designing a controller via the design equations. The 
first test has important practical consequences in terms of 
proper problem setup. It involves direct feed through of 
the disturbance input w in terms of D^: 

F > max ( ar maxt D 1111 D 1 1 1 2 -1 ' CT maxt D 1111 D 1 121 3 ^ * 

The last test involves the characteristics of the solutions 
to the two Riccati equations: 


> 0 
oo 


Y*. > 0 


and 


^max( x » Y ®) < F/ where /i, is an eigenvalue. 

Under the constraints listed and subject to passage of 
the tests an n th order stabilizing controller which achieves 
the inequality | |t| | < T is given by 

D 11 = -Fuji 01 '.!!! (F 2 ! " D 1111 D 111l) D 1112 ~ D 1122 

D 12 and D 21 satisfy 

D , 2 D t 12 = I - D 1121 (r 2 I - D t 1111 D 1111 )‘ 1 D t 1121 and 
F T 21 ^21 = I “ D 1112 (r I - D nil D lin ) d 1112 
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from which 


+ ^12 ) ^12 

C2 — ” 031(^2 ^ 12 ^^ 

B 1 — “H2 + B2D ^2^11 

Ci — F 2 Z + D^D ^21^*2 

A = A + HC + B 2 D _1 12 c i 


and 

z = (I - r' 2 Y 0 xj *’ . 


An n th order controller which achieves the norm bound is 
then given by the realization (A, B 1 , C.,) . There exist 
other controllers which satisfy the norm bound; however, 
they are not necessarily of n th order. The controller 
design equations for this case are omitted. These more 
complex design equations may be found in Glover and Doyle [1J . 


2.0 Transformation to the Standard Form 


In Section 1.0, it was stated that the use of the 
design equations place constraints on the form of D 12 and 
D 21 . This section outlines in detail a method for achieving 
the required form for these matrices. The first step is to 
write them in terms of their singular value decompositions 
as is shown for D 12 : 

°12 = C U 121 U 1 22 H S 12 0] T V T 12 
D 12 = C U 122 U 12lH° S 12 ] T V\ 2 
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D 12 = [U 12 2 u 121 ][0 I]'z'\ 2 v\ 2 . 
Similarly, 

D ?1 = U 2lt S 21 0 3 [ V 211 V 212] 

D ?1 = U 21 [0 I] [V 212 V 2n ] T . 

Now let: 

X, » r 

X '1 = U 21 Jj 21 
u 1 12 = [u 122 <J 121 ] 

^ 1 21 = t V 212 V 211 ] • 

Then 

■>: “ U 'l 2 [0 

0 : , =-• X 21 [0 I ] (v' 21 ) t . 

The equations for z and y are 
z •-= C,x + D ri w + D 12 u 
y — C 2 x ^ D,.w + D ?2 u 
or 

z = c,x + D n w + U' 12 [0 I]X 12 


ORIGINAL PAGE IS 
OF POOR QUALITY 


XVI 


9 



ORIGINAL PAGE IS 
OF POOR QUALITY 

y = C 2 x + X 21 [0 I] (V 1 21 ) T w + D 22 
or 

(U' 12 ) T Z = (U ' 12 ) T C 1 x + (U' 12 ) T D„W + [0 I] T X 12 U 

x_1 2 iY = + [0 I ] (V 1 21 ) T w + X'^D^u. 

Letting 
z' = (U') T z 
u ' = X 12 u 
y- - X'^y 
w 1 - (V' 21 ) t w 
g i ves 
A ' ~ A 

B 'l = S 1 V ' 2 1 

B ' 2 " B 2 X '',2 
C', = (U' 12 )' 

D 'l1 ’ 12^ D 1 ' 21 

d- 12 = [0 I] T 
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C’ 2 = * Va 
D ' 21 = [0 I] 

D ' 22 ~ X 21 D 22^ 12 


The significance of the transformation is that the norms of 
z and w are preserved. This means that designing a 
controller to achieve a particular | |t| |„ for the 
transformed system is equivalent to designing a controller 
to achieve the same goal for the untransformed system, once 
the reverse transformation is applied to the controller. 

i . o Mo difications for Sampled-Data Systems 

More traditional controller design techniques must be 
developed separately for sampled-data and continuous-time 
systems. Fortunately, this is not necessary for H" designs 
due to the fact that the performance criterion has a 
relatively simple frequency domain representation. The 
approach is equivalent to "w-plane" design for single-input, 
single-output systems. The open loop system is assumed to 
have the form 


x(k+l) = Ax(k) 1- B 1 w(k) + B 2 u(k) 
z(k) - qxfk) + D 11 w(k) + D 12 u(k) 
y(k) -• C 2 x(k) + D 21 w(k) + D 22 u(k). 

This discrete time representation can be obtained using 
standard techniques such as those found in Kuo [31 . The 
transfer function matrix is given by 

G(z) = D + C(zl - AJ^B 

where for simplicity B, C, D are appropriate concatenations 
of the open loop system matrices. The bilinear transform 
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z = 1 + w 
1 - w 

is applied to G(z) to obtain 

G ( w) = D + C[ (I+wI) (I-Wl) - A] _1 B 

= D + C(I-Wl) [W(I+A) - ( A-I ) ] _1 B 

= D + C(I-Wl)[Wl - ( I+A) 1 ( A-I ) ] 1 ( I+A) ^B. 

Using the identity 

[l-wl] [wl - ( I+A) 1 ( A-I ) ] 1 = -I + [I - (I+A) - 1 (A-I) ] 

x [wl - (I+A) 1 (A-I) J' 1 

G ( w) = D - C(I+A)'’b + C ( I+A) 1 [Wl - (I+A) 1 (A-I) ] 1 (I+A) "'b 

so that the w-plane state space representation of G(w) is 

D ;> - D - C (I+A) ’b 

C s = 2C ( I+A) 1 

B s = (I+A)' 1 

A s = (I+A) 1 (A-I) . 

The controller is then designed using this representation to 
obtain K(w) , which is represented as a set of continuous 
time state equations with matrices A k , B k , C k , D k . The 
inverse transform 

v/ - z -l 
z + 1 
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is applied to obtain K(z) . State space formulas for this 
transformation are 


A z = -(A k +I) (A k -I) _l 
B z = (I+A z )B k 
C z - C k (I+A k ) 

D z = D k + C k (I+A k ) B k . 

4 . 0 Problem Setup for IMC Controller Design 


The ACES configuration of the Marshall Space Flight 
Center Large Space Structure Ground Test Facility is shown 
in schematic form in Figure 1. The image motion 
compensation (IMC) subsystem is comprised of the line-of- 
sight (LOS) detectors and the IMC pointing gimbals. The 
base excitation table (BET) is the excitation device. 

The natural setup for IMC controller design to minimize 
the effects of BET excitation is outlined here. The 
equations are obtained from a FEM model of the LSS ACES 
ground test facility. The first step in the problem setup 
is to define the z, w, y, and u. Since the actuators and 
sensors are limited to the IMC components, y and u are the x 
and y axis detector signals and pointing gimbal torques, 
respectively : 


u - [IMC X IMC y ] T 
y - [DET X DET y ]\ 

The disturbance vector is most naturally chosen to be the x 
and y axis BET excitation table forces: 


w = [BET* BET y ] . 
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1 . Base Excitation Table 

2. 3 Axis Base Accelerometers 

3. 3 Axis Gimbal System 

4. 3 Axis Base Rate Gyros and Counterweight 

5. 3 Axis Tip Accelerometers 

6. 3 Axis Tip Rate Gyros 

7. Optical Detector 

8. Mirrors 

9. Laser 

1 0. 2 Axis Pointing Gimbal System 

11. LMED System 


Single Structure Control 
Laboratory 


m 



Light Path 


3 Meter Antenna 


Figure 1 . ACES Configuration of the Marshall Space Flight 

Center Large Space Structure Ground Test Facility 
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The first elements of the performance oriented controlled 
variable vector are 


z 1 = [DET X DET y ] T 


since the goal of this system is to reduce the displacement 
of the line-of-sight of the laser beam from its equilibrium 
position on the detector. However, this is not sufficient 
to insure that the conditions of the design equations are 
met. The additional requirement of control input weighting 
is achieved by letting 


2 2 = [IMC X IMC y ] T . 


This satisfies the rank condition on D 12 . The remainder of 
the parameters of the problem can be chosen as follows. 


B 1 is comprised of the appropriate modal gains at the BET 
actuators . 

B 2 is comprised of the appropriate modal gains at the IMC 
actuators . 

C, is comprised of the LOS gains at the detector 
C 2 is also derived from the LOS gains at the detector. 


Also, 


A is in block 2x2 diagonal form and 


D n = [0 0]\ D 12 = [0 I] T 


D 21 0 , D 22 0 . 


22 


The design equations still cannot be applied due to the rank 
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condition on D 21 . This condition is equivalent to requiring 
that the same disturbance enters at two physically separated 
points in the system. The reason for the condition is a 
mathematical technicality. Unfortunately, it places actual 
constraints on the formulation of the problem. 

A possible solution to this rank problem is to define 
D 21 to be a very small constant with respect to the norm of 
the transfer function matrix G 12 at all frequencies of 
interest. 

However, it turns out that this is not the only 
theoretical problem with the above formulation. Another 
difficulty is the minimality condition on the G 12 
realization. The difficulty in the present setup is that 
the requirement is equivalent to the requirement that the 
disturbance (in this case the modes ehich can be excited by 
the BET) must be controllable at the IMC pointing gimbals. 
Unfortunately, this is not the case. In fact, the pointing 
gimbals have significant authority over only four or five 
modes. This is especially troublesome if the technique is 
used without regard to the minimality condition, as the 
design equations will yield a controller without regard to 
the satisfaction of the requirement. In this case, however, 
the controller will not satisfy the norm bound. 

The minimality condition is an artifact of the 
particular procedure used to derive the design equations. 

In the usual factorization approach to H" control, the 
minimality condition is not required since it is possible to 
carry along completely unrelated realizations for each of 
the four transfer functions matrices. It is interesting to 
note that an equivalent requirement for an LQG approach 
would be that the disturbance states be controllable as well 
as observable. 

The consequences of violating the minimality condition 
are illustrated in Figures 2 and 3. Figure 2(a) is the 
magnitude of the open loop frequency response from the x- 
axis BET force input to the x-axis detector. The mode at 
.15 hertz is the AGS hinge point pendulum mode and is 
uncontrollable at the IMC gimbals. The other modes are 
controllable. Figure 2(b) is the magnitude of the closed 
loop frequency response from the x-axis BET force input to 
the x-axis detector. It is apparent that the modes which 
are controllable at the IMC gimbals have been effectively 
suppressed. However, the AGS hinge point pendulum mode 
continues to predominate, although the low frequency 
baseline has been reduced. Figure 3 illustrates the fact 
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that the closed loop norm is not improved. In Figure 3(a) 
the open loop frequency response from the y-axis BET force 
input to the y-axis detector has a maximum of roughly -3 
decibels, as does the closed loop response of Figure 3(b). 
In each case, the values are reliable indicators of the 
infinity norm, since the the x and y axes are only slightly 
coupled. All frequency responses are in transformed 
input/output coordinates, as discussed in Section 

As these problems were uncovered in the attempt to 
obtain an IMC controller design, it was decided to 
investigate the properties of the design equations using an 
extremely simple model. The next section documents the 
findings of this investigation. 
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Figure 3(a) . Open Loop Frequency Response Magnitude 
(BET y to LOS y ) 



Figure 3(b). Closed Loop Frequency Response Magnitude 
(BET y to LOS y ) 
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5.0 


SIMPLE TRACKING PROBLEM 


The usual way to introduce design flexibility in the 
factorization approach to H* optimal control is via 
frequency dependent weightings. To see the effect of 
frequency dependent weighting, the usual parametrization of 
the closed loop transfer function is useful. For a stable 
plant, the closed loop transfer function can be written as 


T, - T 2 QT 3 
where 



Any stable transfer function Q generates a stable closed 
loop transfer function and a controller which achieves that 
transfer function. In fact, every stable closed loop 
transfer function is generated by some stable Q. 

The usual requirement for a solution to exist is that 
T^ and T 3 have constant rank on the extended jw axis. No 
minimality condition is required. Weighting is introduced 
by solving the modified problem of minimizing the infinity 
norm of 


t, - t 2 w 2 qw 3 t 3 


where W 2 and W 3 are chosen to satisfy the constant rank 
condition and to define frequency ranges over which 
optimality is emphasized. The major difference in the 
general H problem and the problem solved by the new design 
equations in question is the presence of the two minimality 
conditions. The effects of the constraints are most easily 
seen by examining the block diagram of Figure 4. 
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Figure 4. Block Diagram Used for Examining the Effects 
of Frequency Dependent Weighting. 

The minimality constraints imply the following: 


The system realization is simultaneously 
completely controllable by u and observable 

by z . 


The system realization is simultaneously 
completely controllable by w and observable by 
Y- 


The first can be satisfied by including the signal y in 
the definition of z. A consequence of this is that T, is 
also weighted by W 3 . The second condition can only be 
satisfied if W 2 has zero dynamical order and is a 
fundamental limitation of the new algorithm. It should be 
pointed out that if W 2 is derivative in nature, an 
appropriate realization can be obtained, but it will 
generally increase the order of the plant, and hence the 
controller. Also, the appropriate modifications to account 
for nonzero D 22 must be used. The actual details of 
combining the weightings into an appropriate realization 
remain to be worked out. 

Although no results are available using the weighting 
scheme of Figure 4, the alternate weighting scheme of Figure 
l > was used with an extremely simple plant to illustrate what 
is possible when all of the constraints of the new algorithm 

are satisfied. 
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Figure 5. Alternate Weighting Scheme Used for a Simple 
Tracking Problem 

The plant transfer function is 

G (s) = 1/ (s + 1) . 

The weights are given by 

W, == . 99 ( s/ 10 + 1) 

W-, = . 001 . 

The z vector is defined by 


z i = Yg " Yw2 
7*2 - 10 5 U . 


It is interesting to note that although the plant is scalar, 
the two-dimensional nature of the z vector means that H* 
control is inherently a multivariable problem. The z 2 
element is chosen small enough to simplify the process of 
obtaining an approximate solution. Another important point 
to note is that the .99 multiplier in W, is necessary to 
achieve a high gain controller. 


Figure 7 is the closed loop frequency response when the 
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controller is implemented in the block diagram of Figure 6. 
The important point to be made here is that an H* approach 
can be used to design a simple tracking system with 
specified closed loop bandwidth (by choosing the break 
frequency of W,,) and specified steady-state error constants 
(in this case less than .01 error to a unit step input). 


CONCLUSIONS AND RECOMMENDATIONS 


The design equations and appropriate constraints and 
assumptions for a new simplified algorithm for designing H” 
optimal controllers have been reviewed. A transformation 
required to satisfy an important constraint has been 
cl rived . The use of the design equations with sampled-data 
systems is outlined in detail and the required state space 
transformations are summarized. 

The requirement of using frequency-dependent weights is 
illustrated via two examples. One of the examples is a 
simplified but realistic design for a subsystem of a large 
space structure ground test verification facility. The 
other is used to illustrate the efficacy of a particular 
weighting scheme. 

further work is indicated in the area of developing 
stategies for finding weights to achieve particular design 
goals. It is also suggested that possible tradeoffs between 
performance and robustness be investigated by studying 
various problem formulations. 
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